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A single-photon Fock state has been generated by means of conditional preparation from a two- 
photon state emitted in the process of spontaneous parametric down- conversion. A recently de- 
veloped high-frequency homodyne tomography technique has been used to completely characterize 
the Fock state by means of a pulse-to-pulse analysis of the detectors' difference photocurrent. The 
density matrix elements of the generated state have been retrieved with a final detection efficiency of 
about 57%. A comparison has been performed between the phase-averaged tomographic reconstruc- 
tions of the Wigner function as obtained from the measured density-matrix elements and from a 
direct Abel transform of the homodyne data. The ability of our system to work at the full repetition 
rate of the pulsed laser (82 MHz) substantially simplifies the detection scheme, allowing for more 
"exotic" quantum states to be generated and analyzed. 

PACS numbers: 03.65.Wj, 03.65.Ud, 42.50.Dv 
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INTRODUCTION 



Homodyne detection is a widely used method to obtain 
tomographic reconstructions of classical and quantum 
optical fields in the cw and pulsed regimes. It consists in 
mixing the unknown signal field with a strong local oscil- 
lator (LO) field on a beam splitter whose outputs are then 
detected by proportional photodetectors 'l'|. The differ- 
ence in the photocurrents produced by the two detectors 
is proportional to the electric field quadrature selected 
by varying the relative phase between the LO and the 
signal field. From a collection of field quadratures it is 
then possible to reconstruct the quantum state by means 
of tomographic techniques. Vacuum, thermal, coherent 
and squeezed states have been successfully analyzed, in 
the frequency-domain, by means of a spectral analysis of 
the difference photocurrent y, |3|, and in the time do- 
main, with the measurement of integrated photocurrents 
on isolated pulses, in order to achieve a full reconstruc- 
tion of the unknown field states jj, y, 1^1 . 

Recently, a pulsed homodyne technique 7] has been 
successfully applied to the first complete quantum state 
reconstruction of single-photon Fock states of the elec- 
tromagnetic field 1^, i.e. states containing exactly one 
quantum of field excitation. Such states are an impor- 
tant workbench for the test and the analysis of nonclassi- 
cal behavior and represent a fundamental building block 
for the implementation of efficient linear quantum com- 
putation. Time-domain homodyne detection schemes as 
those used in m, Iq are however quite demanding from 
the experimental point of view, since they require very 
low electronic noise and high subtraction efficiencies over 
very large frequency bandwidths (from dc to a few times 
the repetition rate of the laser) . Mostly due to these tech- 
nical limitations, all the existing applications of the lat- 
ter technique have been limited to pulse repetition rates 



well below the megahertz. The measurements reported 
in |3, y indeed made use of a pulse picker to lower the 
laser repetition frequency, and the low state production 
rates implied very long acquisition times in order to store 
enough data for a single accurate tomographic recon- 
struction. More recent works by the same group [sl, ll(l | 
have adopted a different technique, where the local os- 
cillator pulses are sent "on demand" to the homodyne 
detector only when needed; this has contributed to sig- 
nificantly increase their acquisition rates, although at the 
price of a more complicated setup and lower efficiencies. 



Our group has recently developed an homodyne set- 
up which is able to perform high-frequency (with repe- 
tition rates up to about 100 MHz) time-domain analy- 
sis of quantum states of light [ll|. Its capabilities were 
first demonstrated by performing a tomographic recon- 
struction of the Wigner function of weak coherent states 
and by showing the possibility of a gated-mode acqui- 
sition for conditional measurements. Here we present 
the first complete and efficient tomographic reconstruc- 
tion of a single-photon Fock state based on such a time- 
domain, pulsed optical homodyne apparatus operating 
at the high repetition frequencies characteristic of com- 
monly used mode-locked laser systems. By allowing us 
to perform ultrafast pulse-selective measurements, this 
apparatus exhibits two fundamental advantages over al- 
ternative schemes working at lower repetition rates: on 
one side, thanks to the much shorter acquisition times, 
the stability requirements of the overall experimental sys- 
tem are loosened; on the other, if long-term stability 
of the setup is available, states exhibiting stronger non- 
classical features (usually characterized by smaller gen- 
eration probabilities) can be efficiently analyzed. 



GENERATION OF SINGLE-PHOTON FOCK 

STATES 

A single-photon Fock state can be generated by means 
of a conditional state preparation on the two-photon 
wavefunction emitted in the process of spontaneo us p ara- 
metric down-conversion (SPDC) (see Refs. [ij, [iSJ and 
references therein). Such a state is written as 



I*) = \0)J0)^- d^k,duJsd^ktdu;t^j{k„LU„kt,LUt) 

X \ks,UJs)s\^t,^t)t, (1) 

where the function ip(\i.s,LUs,kt,LL!t) describes the spatial 
and spectral properties of the two-photon state {bipho- 
ton) and the subscripts s and t correspond to the signal 
and trigger (idler) photons. 

In order to non-locally select a pure state on the sig- 
nal channel, trigger photons must undergo narrow spatial 
and frequency filte ring before being detected by a single 
photon counter |l2tll4l |. The conditional non-local prepa- 
ration of the signal state is thus described by the density 
operator ps given by the partial trace of the product be- 
tween the SPDC-state density operator and the measure- 
ment density operator pt- 



where 



Pt 



p, = Tr{p,|*)(*|}t (2) 



d^ktduJtT{\<it,ujt) \kt,u}t)^ (kt,wt|j (3) 



and ^(k, lo) is the transmission function of the spatial 
and spectral filter placed along the trigger photon path. 
Note that the nonlocally-prepared signal state will only 
approach a pure state if the filter transmission func- 
tion T(k, w) is much narrower than the momentum and 
spectral widths of the pump beam generating the SPDC 
pair. Accordingly, one can define a purity parameter as 
P = Tr(/5^) which, in the spectral domain, is expressed 
as 



Pf = 



1 



-)l-l 



(4) 



where Gaussian profiles with widths Up and a f are as- 
sumed for the pump power spectrum and for the trigger 
filter transmission, respectively. It is easy to see how 
the spectral purity parameter Pt approaches unity when 
Uf ^ dp. In such a condition the strong spectral filter- 
ing performs a non-local selection of a pure signal state 
whose properties are defined by the pump as shown in 
[l.! Illl. Similarly, the spatial purity parameter is found 
to be 1121: 
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where Ki and Kp are the beam widths in the momentum 
space for the idler spatial filter and for the pump, respec- 
tively. 



EXPERIMENTAL SETUP 

A mode-locked Ti:sapphire laser, emitting 1-2 ps long 
pulses at 786 nm with a repetition rate of 82 MHz is 
used as the primary source for the experiment, schemati- 
cally drawn in FigQ] The laser pulses are first frequency 
doubled in a 13-mm long LBO crystal which produces 
the pump pulses centered at 393 nm, with a mean power 
of about 100 mW. The UV pump pulses, after a spa- 
tial mode cleaner, are then slightly focused, with a beam 
waist Wp — 220 /im, and the waist position is carefully 
located inside a 3-mm thick, Type I BBO crystal cut for 
coUinear degenerate SPDC. The crystal is adjusted in or- 
der to obtain a signal-trigger cone beam with an angle of 
^ 3°. Signal and trigger photon pairs are then selected 
by means of irises placed at 70 cm from the output face 
of the crystal. 
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FIG. 1: Experimental apparatus: BBO downconverter crys- 
tal, PZT piezoelectric transducer, PBS polarizing beamsplit- 
ter, A/2 half-wave plate, HD homodyne detector. 

The trigger beam is directed to the state preparation 
channel, where it is spectrally filtered by means of a pair 
of etalon interference filters and is then coupled into a 
single-mode fiber connected to a single photon counting 
module with an active-area diameter of 200 /im (Perkin- 
Elmer SPCM AQR-14). The typical trigger count rate 
after spectral and spatial filtering is of about 500 s""'^, 
with a dark count rate of the order of 10 s^^ when both 
are gated with the laser pulses. 

The combination of the two etalon cavities yields a fi- 
nal spectral width for the trigger channel of about 50 
GHz. By comparing this with a measured pump width 
of 430 GHz, we estimate a spectral purity parameter of 
Pt — 0.98. However, due to the second harmonic gener- 
ation process, the LO spectral width remains ^/2 times 
narrower than that of the selected single-photon state. 



The degree of mode- matching between the LO and the 
single-photon state cannot be improved by filtering the 
LO beam and might be increased only by inserting a 
spectral filter on the UV pump. However, this solution 
would also drastically lower the count rate on the trigger 
detector and has not been implemented yet. 

Filtering in the spatial domain is achieved by means of 
the single-mode fiber and our experimental parameters 
lead to a spatial purity parameter P^ = 0.86, which again 
means that the conditionally prepared single-photon 
state mode is essentially determined by the pump char- 
acteristics. The LO beam is spatially mode-matched to 
the fiber-selected signal photon mode by the insertion of 
appropriate lens combinations along its path. In order 
to finely adjust the alignment and the synchronization 
between the signal and LO pulses, we use the stimulated 
beam produced by injecting a seed pulse into the crystal. 
Under appropriate conditions [l^, the amplified seeded 
beam can spatially well simulate the selected signal beam 
and can be used for alignment purposes. 

The LO pulses are obtained by splitting a small portion 
of the laser emission, and their polarization is rotated 
by a half-wave plate in order to overlap them with the 
conditionally-prepared single-photon signal beam onto a 
polarizing beam splitter (PBS). One of the steering mir- 
rors is mounted on a piezoelectric transducer in order to 
vary and control the LO phase. A 300-mm focal-length 
lens, an additional A/2 wave plate, and a second polariz- 
ing beam splitter are then used to accurately mix the two 
field modes and focus them onto the two photodiodes of 
the homodyne detector. 



HIGH-FREQUENCY BALANCED HOMODYNE 
DETECTION 



Two photodiodes (Hamamatsu S3883, with active area 
1.7 mm^) are connected to the positive and negative in- 
puts of an operational amplifier providing a gain of about 
10. The bias voltage of the two photodiodes is optimized 
in order to provide the maximum extinction in the dif- 
ference signal when the optical power impinging on the 
two is well balanced. The amplifier output is AC-coupled 
(cut-off frequency of 1 MHz at -3 dB) to a second am- 
plifier with gain G =100 and bandwidth extending from 
dc to 300 MHz, before being sent to a digital scope with 
an analog bandwidth of 1 GHz, a sampling rate of 10 
Gsamples/s and a vertical resolution of 8 bits. 

In order to decrease the effect of the dark counts in the 
single-photon trigger detector, a strict coincidence with 
the signal coming from the laser pulse train is used as 
the trigger for the acquisition of the homodyne signal. 
Although this slightly reduces the trigger count rate, it 
is effective in increasing the ratio of "true" to "false" 
trigger events to more than 99%. 

A single temporal frame composed of 250 acquisition 



points is stored in the oscilloscope memory whenever a 
trigger event occurs. Each of these frames contains two 
consecutive LO pulses where only the first one is coin- 
cident with the detection of a trigger photon and con- 
tains the "single-photon information", while the second 
one can be used for the measurement of the reference 
vacuum state. About 5000 frames can be stored in a 
sequence at a maximum rate of 160,000 frames per sec- 
ond. With an average rate of trigger counts of about 300 
s~^, this means that each sequence is acquired in about 
15-20 s. Each sequence of frames is then transferred to 
a personal computer where the areas of the pulses are 
measured and their statistic distributions are analyzed 
in real time. The acquisition of two nearby LO pulses in 
each frame allows the simultaneous analysis of the "sig- 
nal" pulses, namely those that have interfered with the 
signal photons, and of the "dark" pulses, so that both 
the Fock-state and the LO shot-noise distributions are 
immediately available from a single sequence, indepen- 
dent from low-frequency fluctuations in the system. The 
acquisition and analysis of a sequence of 5000 signal and 
dark pulses yields an experimental Wigner marginal dis- 
tribution in just about 30 s. 

It is interesting to compare these acquisition times with 
those reported in |3| by the Constance group for the first 
tomographic reconstruction of the single-photon Fock 
state. Due to the limited bandwidth of their homodyne 
detector (1 MHz), the overall pulse rate had to be low- 
ered to about 800 kHz by means of a pulse picker at the 
exit of the laser. In that case about 12,000 experimental 
data points were acquired in a 14-hour experimental run 
in order to obtain an accurate phase-averaged marginal 
distribution of the quantum state. Our setup now al- 
lows the same kind of acquisition to be performed in just 
about one minute, with a gain of almost three orders of 
magnitude in the measurement time. It is clear that such 
a speed-up can help to loosen some of the constraints on 
the overall stability of the generation system and, in par- 
ticular, on the interferometric stability of the beam paths 
which is needed when phase-dependent quantum states 
are to be analyzed. On the other hand, by greatly in- 
creasing the overall number of available LO pulses in a 
reasonable experimental time interval, our scheme will al- 
low us to investigate much rarer and more exotic events, 
such as those involving higher photon numbers |llL Il6| . 

Note that recent experiments from the same group 
[Sl lUl have seen the adoption of a clever solution to 
the limited acquisition rate caused by the low detector 
bandwidth. Instead of reducing the whole experimental 
repetition rate, the pulse picker has been placed only in 
the LO path and activated "on demand" upon a trig- 
ger detection event. This has allowed to keep reason- 
able pair production rates and to obtain a final acqui- 
sition rate similar to the one presented here. However, 
the introduction of the pulse picker, besides significantly 
complicating the experimental setup, has the important 



drawback that the trigger delay of the shutter in the LO 
path has to be compensated by a long (about 15 m) op- 
tical delay line on the signal photon path in order for 
the two to reach the mixing beam-splitter at the same 
time. This clearly introduces additional losses and may 
degrade both the mode quality and the accurate phase 
stability of the fragile quantum state. 

The performances of our detector were verified with a 
preliminary measurement of the laser shot-noise by ana- 
lyzing the variances of the "dark" pulse areas as a func- 
tion of the LO power. A linear behavior is preserved 
over a wide range of LO powers, extending up to 9 mW, 
and with a signal to electronic- noise ratio of about 12 dB 
when the device is operated at the optimum LO power 
of 7 mW. 



QUANTUM STATE RECONSTRUCTION 

Balanced homodyne detection allows the measurement 
of the signal electric field quadratures xg as a function 
of the relative phase imposed between the LO and the 
signal. By performing a series of homodyne measure- 
ments on equally-prepared states it is possible to obtain 
the probability distributions p{x, 9) of the quadrature op- 
erator Xg that correspond to the marginals of the Wigner 
quasi-probability distribution Wix^y) |l7l |: 



p{x,e) 



+ 00 



W{x COS 6 — y sin 6, x sin 6 + y cos 6)dy. 



(6) 
In principle, given a sufficient number of quadrature 
distributions for different values of 0, one is able to re- 
construct the Wigner function by inverting Eq. ©. In 
the case of Fock states, some simplification in the recon- 
struction procedures can be accomplished thanks to the 
phase-invariant nature of the Wigner function 8] . Only 
one phase-averaged marginal distribution is required for 
the reconstruction in this case and, from the experimen- 
tal point of view, this also means that the stabilization of 
the LO phase is not needed. Furthermore, the Abel trans- 
form, relying on the cylindrical symmetry of the problem, 
can be used instead of the Radon transform in order to 
reconstruct a section of the Wigner function W{r) [ig, 
as given by 



W{r) 



-t-oc 



dpa.-v{x) dx 
dx 



V^ 



(7) 



where Pa.-v{x) is the phase-averaged marginal distribu- 
tion. The numerical implementation of this method is 
however extremely sensitive to noise and a large num- 
ber of data is required to obtain a reliable result. Since 
this reconstruction procedure involves the application of 
numerically-unreliable methods, the group of D'Ariano 
has developed an alternative precise technique that per- 




FIG. 2: Vacuum (open circles) and signal (filled circles) phase- 
averaged quadrature field distributions. The solid fines are 
fitted curves. The quadrature is normalized to the vacuum. 



mits the reconstruction of the density matrix elements 
directly from the homodyne data |l9l |. 

The distributions of the acquired pulse areas for the 
vacuum and the single-photon state are shown in Fig. |21 
where each histogram is obtained from the analysis of 
about 200,000 acquired pulses. Such distributions rep- 
resent the marginal integrals of the corresponding field 
state. The non-unity detection efficiency of the appa- 
ratus prevents us from observing the real single-photon 
Wigner function, and what we get instead is its convolu- 
tion with the vacuum one. The convolution result is the 
well known s-parametrized quasi-probability distribution 
with the s parameter scaled by the detection efficiency 
T] [lalSdl- The expression for the corresponding marginal 
distribution is: 



p{x;r)) 



[1-77(1 



4a;2)le-2-^ 



(8) 



From a fit of the experimental distributions to the cor- 
responding theoretical curves we obtain an overall detec- 
tion efficiency of ?? = 0.574 ± 0.002. 

The expected overall efficiency rj is given by the con- 
tribution of several terms: 



V = VhdridcVn 



(9) 



where ryhd = 0.90 is the efficiency of the homodyne de- 
tector including the PBS losses, while ride = 0.99 is 
the efficiency connected to the dark counts of the trig- 
ger detector. The mode-matching efficiency is given by 
Vmm = ilexpVPsPt, whcrc \/ PsPt depends on the tight- 
ness of the spatial and spectral filtering on the trigger 
channel and constitutes the upper limit for the mode 
matching efficiency achievable between a pure state (LO) 
and the non-locally selected single photon state. The 
term rycxp is connected with the non-ideal experimental 
conditions and can be estimated from the visibility of 
the interference fringes formed between the LO and the 
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FIG. 3: Reconstructed diagonal elements of the density ma- 
trix for the single-photon state and the vacuum (inset). 



stimulated radiation used for alignment ilj|. We obtain 
77exp — 0.7 and the overall efRciency is thus expected to 
be 77 ~ 0.6, in good agreement with the value retrieved 
from the fit of the marginal data to the theoretical dis- 
tributions. 

In Fig. 131 the reconstructed diagonal elements pn,n of 
the density matrix are reported as obtained by averaging 
the so-called pattern functions over our data: 



Pn,n = \n\p\n) 



4-CX2 



Pav(a;)/„„(a;)dx. 



(10) 



Here /„„(a;) are the factorized pattern functions that 
can be numerically implemented following the numer- 
ical recipes given in Ref. [1^ |21|. Figure El clearly 
shows the non-poissonian photon distribution of the re- 
constructed state, which is a statistical mixture of the 
single-photon Fock state with the vacuum state p = 
Pofi |0) (0| + pi.i |1) (1| with po,o = 0.426 ± 0.001 and 
pi^i = 0.572±0.002 in good agreement with the efficiency 
value found by fitting the mar gina l distributions. 

As discussed by D'Ariano in |2l|,l22?|, the Wigner func- 
tion can be obtained from the reconstructed density- 
matrix elements in a dimensionally-truncated Hilbert 
space as: 



W{x,y)=ReY. 
and 



M M-d 

id arctan(y/a;) 



d=0 



^ A{n,d,x,y)pn.n+d 

(11) 



n=0 



A{n,d,x,y) = ^^^[2-SM]\2{x + ty)f 



where Lf^ix) are Laguerre polynomials, S is the Kro- 
necker delta and M is the maximal quantum number of 




FIG. 4: Wigner function of the single-photon Fock state as ob- 
tained from the reconstructed density-matrix elements. The 
negativity of the distribution, a clear proof of the non-classical 
character of the state, is evident around the origin. 



the reconstructed density matrix. In Fig0]we show the 
Wigner function obtained from the ffi'st 10 reconstructed 
diagonal elements of the density matrix and by setting 
all the off-diagonal terms to zero thanks to the phase- 
invariance of the state. The reconstructed Wigner func- 
tion assumes negative values at the origin, showing the 
non-classical features of the single-photon state. We find 
W^moas(0,0) — —0.0995, which is quite close to the theo- 
retical value of the Wigner function W{0, 0; 77) — —0.0942 
obtained by using Eq. Hll|l with the mixed state defined 

by 



/5=(l-,7)|0>(0|+rHl)(l| 



(13) 



with the efficiency value (77 = 0.574) retrieved from the 
fit of the phase-averaged marginal distributions. 

A comparison between the different methods of Wigner 
function reconstruction is presented in Fig. [S] where the 
W{x,0) section is plotted. Both the proffle obtained by 
a direct implementation of the Abel transform on the ac- 
quired homodyne histograms using the basic algorithm of 
Nestor and Olsen |2^ and a section of the Wigner func- 
tion as reported in Fig||] are shown, together with the 
curve retrieved using the mixed state of Eq. H13() . It is 
evident that the direct reconstruction method based on 
the Abel transform is extremely sensitive to the noise 
present in the data and fails to accurately reproduce 
the expected vacuum-convoluted single-photon Wigner 
function. On the other hand, the reconstruction via the 
density-matrix elements is much smoother and extremely 
close to the expected proffie. 

Recently, alternative reconstruction techniques based 
on maximum-likelihood estimation methods have been 
proposed |2J, [23 and used to analyze experimental 




FIG. 5: W{x,0) section of the reconstructed Wigner func- 
tion. Filled circles: Abel transform; empty squares: section 
of Fig|l] solid line: section corresponding to the mixed state 

of Eq. rra . 



data !lC| . Such methods can reconstruct the density ma- 
trix elements and the Wigner function with lower require- 
ments on the number of data samples without a substan- 
tial increase in the statistical errors. The implementation 
of these methods might contribute to further decrease 
the acquisition time of our experiment whenever this be- 
comes a limiting factor, i.e. when more exotic quantum 
states are analyzed. 



CONCLUSIONS 



We have applied a recently developed technique of 
time-domain, high-frequency homodyne detection to per- 
form quantum tomographic reconstructions of the single- 
photon Fock state. The scheme for the generation of 
the quantum field is based on remote state preparation 
from an entangled two-photon wavefunction emitted in 
the process of pulsed parametric down-conversion. By 
analyzing pulsed homodyne signals at the full repetition 
rate of the mode-locked source laser, this scheme allows 
us to perform full tomographic reconstructions in rela- 
tively short times. Furthermore, by allowing us to avoid 
the introduction of additional mechanisms of pulse selec- 
tion, our scheme has the potential to achieve more com- 
pact realizations (which may be an issue when an accu- 
rate control of the relative phase in needed) and higher 
detection efficiencies in the analysis of quantum states. 
Finally, we have shown that the methods of Wigner func- 
tion retrieval based on the direct analysis of the homo- 
dyne data are too sensitive to experimental noise and 
that much more accurate results can be obtained from 
the reconstructed density-matrix elements. 
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